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“II n’est pas necessaire d’esperer pour entreprendre, 
ni de reussir pour perseverer. ” 
(Begin, even without hope; Proceed, even without success.) 

William of Orange (“William the Silent”), murdered in Delft in 1584. 

Abstract. Difference calculus compatible with polynomials (i.e., such that the 
divided difference operator of first order applied to any polynomial must yield a polyno¬ 
mial of lower degree) can only be made on special lattices well known in contemporary 
q— calculus. Orthogonal polynomials satisfying difference relations on such lattices are 
presented. In particular, lattices which are dense on intervals (|g| = 1) are considered. 

1. Introduction. 

Many works have been devoted to orthogonal polynomials satisfying remarkable 
differential or difference relations. 

For instance, the classical orthogonal polynomials are characterized by the exis¬ 
tence of a differential relation of the form 

W(x)p' n (x) = UJ n (x)p n (x) + d n p n -i(x), (l) 

where W is a fixed polynomial of degree < 2, u> n is a polynomial of degree < l,and 
where d n is constant (in x) (cf. [1], p.8). 

This differential relation interacts most efficiently with the recurrence relation 

a n+1 p n+1 (x ) = (x- b n )p n (x) + a n p n _!(x) (2) 

in the production of various useful identities. 
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If we accept higher degrees polynomials in (1): 


W(x)p' n (x) = fi n (x)p n (x) - a n O n (x)p n - 1 (x), (3) 

with W, Q n and 0 n polynomials of degrees < s + 2, s + 1 and s > 0, we get the semi- 
classical class, studied by Laguerre (the notations of (3) are [almost] Laguerre’s ones 
[19]), Hendriksen and van Rossum [10], and Maroni [22] who coined this name. See 
also [4,5] for determination of the relevant measure. 

If we try to extend (3) to a difference operator of first order, we expect to see 
the derivative p' n (x) replaced by some combination of p n (y(s )) and p n (y(s + 1)), where 
y(s) and y(s + 1) are two consecutive points on a lattice associated to the difference 
operator. 

We will explore here the extension of the semi-classical property (3) to the re¬ 
markable difference operators and the corresponding nonuniform lattices studied by 
many people recently [2,6,12,16,17,21, 24,25,29,30]. In particular, those lattices which 
happen to fill densely an interval will be considered with special care. 


2. The difference operator and the related lattices. 


We consider here a first-order difference operator involving the values of a function 
at two points. For each x, let (fi(x) and P 2 {x) be these still unknown points. The first- 
order divided difference operator at x is 


(Pf){x) 


f(<P 2pg)) - f(<P l(gQ) 

<P2(x) ~ <Pl(x) 


(4) 


If we impose the condition that Vf is a polynomial of degree n — 1 if / has degree n, 
then <fi(x) and <f 2 (x) must be the two roots in y of a quadratic equation 


Ay 2 + 2 Bxy + Cx 2 + 2 Dy + 2 Ex + F = 0. 


(5) 


(see [6,12, 21]). 

Indeed, applying V to f(x) = x 2 and f(x) = x 3 readily yields that + p >2 and 
<pf ~h < ~pi l ~P 2 + P 2 must be polynomials of degrees 1 and 2, which implies (5). Conversely, 
if (5) holds, any symmetric polynomial in cpi and </?2 is a polynomial in x. 

Let us figure the conic (5) and one of its parametric representations {x(s),y(s)} 
such that y(s) and y(s + 1) appear naturally as the two ordinates associated to the 
abscissa x = x(s): one starts from some point {x\ = x(si),yi = y(si)} on the conic, 
and one looks for the points {xi- = x(s 1 + k — 1), y^ = y{s\ + k — 1)}, k — 1,2,... 
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Figure 1. x— and y— lattices. 


To achieve this, let us consider the familiar parametric representations in the two 
following cases: 

1. The conic (5) has a center B 2 — AC ^ 0. With the center coordinates x c = 
(AE - BD)/(B 2 - AC) and y c = (CD - BE)/(B 2 - AC), one has A(y - y c ) 2 + 
2 B(x — x c )(y — y c ) + C(x — x c ) 2 + F = 0, with F = F — Ay 2 c — 2Bx c y c — Cx 2 = 
F + Dy c + Ex c = F + (CD 2 - 2 BDE + AE 2 )/(B 2 - AC), 

X = x(s) = x c + C \fA(q s + q ~ s ) , y = y(s) = y c + (Vc(q s ~ 1/2 + q ~ s+l/2 ) , (6) 

is a valid parametric representation of (5), (if AC ^ 0), where 
(2 = F/(4(B 2 - AC)) and 


qin + q -l/2 = _ 


\/ r AC 


, 4 B 2 

^■’9 + 9 - =Jc~ Z 


Indeed, x(s) = x(—s) is kept by the transformation s <-> — s in (6) but y becomes 
y(s+ 1). 

If AC = 0, the scheme of Fig. 1 and Fig. 2 does not work: horizontal and/or 
vertical lines do not meet the conic in two points any more. 
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Figure 2. Other kinds of lattices. 

The generic (also rightly called hyperbolic) case \q\ ^ 1 gives a hyperbola (Fig. 2(a)). 
Remark that the asymptotes are given by (x and y large) y (C/A) l / 2 q ±l / 2 x: 
q is the ratio of the slopes of the asymptotes of the conic. 

If F = 0, one finds x — x c = X\fAq s , y — y c = X\fCq s±1 ^ 2 , the “old” q— lattice 
(Fig. 2(b)). 

2. The conic (5) has no center, B 2 — AC = 0. Then, 
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found directly to satisfy (5) in the parabolic case B 2 = AC, or by taking the 
limit of (6) when q —> 1: let q = exp(e), then, q s = 1 + se + s 2 e 2 /2 + • • 

B = -VAC{q 1/2 + q~ 1/2 )/ 2 = -\fAC{ 1 + £ 2 /8 +•••), etc. (Fig. 2(c)). 
if dVc + eVa = o, one redefines s through the translation 


Sold 


_lim ± 

dVc+eVa ^o 


Snew 


+ 


V ( Q 2 _ af)c j 

2(F>\/C + ^\/A) J 


? 


to get 


x = —2i 


'D 2 -AF 


E 


AC 


■s ,y = 


v/IC 


- 2 


\[EB^AF 

I 


(a ±1/2) , 


the simplest lattice (this one is uniform (Fig. 2(d))). 

These points form one of the special non uniform lattices ( snul ) I-VI of [24,25]. If the 
quadratic equation (5) describes an ellipse, as in the figure, this lattice fills densely an 
interval, unless it is finite (periodic: q N = 1). 


3. Semi-classical orthogonal polynomials on snul s. 

Semi-classical snul orthogonal polynomials may be defined through a V —difference 
equation of the form 


W(x)(VS)(x) = 2V{x){MS)(x) + U{x) (8) 

for the Stieltjes function 

p OO 

S(x) = / (x — f)~ 1 d/i(f) = y~] Hk/x k+1 

J SUPP. /I Q 

(Stieltjes transform of the orthogonality measure d/i) where W, V and U are polyno¬ 
mials and A4 is the arithmetic mean operator: 

( Mf)(x ) = (f(<p i(x)) + f(<p 2 (x)))/2. 

If id has a jump (pole of S) at some y n , it must have jumps at the other points 
y n ± i, y n ± 2 i- ■ ■ of the corresponding lattice. 

It is possible to recover a difference relation extending (3) to the present difference 
operators, as attempted and min of meer (less or more) achieved in the paper [21], see 
here an attempt to achieve the converse: 

Theorem. Orthogonal polynomials are V—semi-classical, i.e., their Stieltjes function 
satisfies an equation (8) with polynomials W, V, and U, if each p n satisfies a linear 
first-order difference relation connecting p n with p n -\ 

W n (x)(Vp n )(x) = fi n (x)(Mp n )(x) - a n O n (x)(Mp n -i)(x), (9) 
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where W n , Q n , and 0 n are polynomials of fixed (independent of n) degrees. 

Indeed, let y\ and y 2 be the ordinates corresponding to x — aq, then (9) is a linear 
relation involving p n (yi), Pn(y 2 ), p n -i{yi), and p n - 1 (y 2 ): 

A n p n (yi ) + B n p n (y 2 ) + CnPn -1 (j/i) + A»Pn- 1 ( 2 / 2 ) = 0, (10) 

where A n , B n , C n , and are rational functions of fixed degrees of x = x\ and y\ (y 2 
may be replaced by —2 (Bx + D)/A — y\ in A n = — W n (x)/(y 2 — y\) — fl n (x)/2 etc.). 
We come to n + 1: 

A n +ip n +i(yi) + B n+ ip n+ 1 ( 2 / 2 ) + C n +ip n (yi) + D n+ ip n (y 2 ) = 0, 
and use the recurrence relation (2): 

[(yi — bn)A n+ 1 + a n+ iC n+ i] p n (yi) + [(y 2 — b n )B n+ 1 + a n+ iD n+ i] p n (y 2 ) — 

A n -\-\P n — l {x/\ ) dji-Hji+lPn-1 (?/2 ) — 0, (ff) 

which is also a linear relation involving p n (yi), p n ( 2 / 2)5 Pn-i(yi), and p n - 1 ( 2 / 2 )! 

1. Either (10) and (11) are dependent for each n, then A n = K n [(yi — b n )A n+ i + 
a n +iC n+ i] and C n = K n a n A n+ 1 , i.e., the recurrence relation A n = K n [(y 1 — 
b n )A n+ 1 — a^ +1 ET n +i^-n+ 2 ] which looks somewhat like (2): actually, 

K 1 ... ET n -i®i... a n -\A n would be a solution of (2), and this is incompatible with 
the requirement that A n keeps a finite degree. So, 

2. or (10) and (11) are independent for some n. Then, it is possible to extract p n (y 1 ) 
and p n {y 2 ) in terms of p n _ 1 ( 2 / 1 ) and p n - 1 ( 2 / 2 ) ■ 

Pn(yi) = X n p n _ 1 (yi) + Y n p n -i(y 2 ) , Pn{y 2 ) = ^nPn-i(yi) + U n p n - 1 (y 2 ), 

where X n , Y n , Z n , and U n are again rational functions of fixed degrees. We get 
the relation for n + 1 using (2): 


Pn+l(yi) 

1 

yi - b n 

0 

Pniyi) 


Pn- l(2/l) 

Pn+l(y 2 )_ 

1 

0 

y2 — bn 

_Pn(y 2 ). 

®n+1 

.Pn-1(2/2)_ 


/ 1 

2/i b n 

0 


"X n 

w" 

_1 1 

Pn (//l) 

«n+l 

0 

2/2 — b n _ 

®n+l 

_ 

u n _ 

1 

_Pn(2/2). 


showing that we can construct similar relations for n + 1, n + 2, etc. as far as x is 
not a zero of the determinants S n = X n U n — Z n Y n , etc.: 


Pn+1 (2/1) 


Wi+1 

Wi+1 

Pn (2/l) 

Pn+1 ( 2 / 2 ) 


■^Vi+l 

U n -\-\ 

Pn ( 2 / 2 ) 
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®x 


whence the recurrence relations for the X’s, Y’s, Z' s, and C/’s: 

l-^n+l — V 1 b n Q*nUn/i 
®n+lC/n+l ' V2 — b n dnX n /8 n , 

®n+l^n+1 8 n , 

®n+l^n+l — ® nZ n /8 n j 

fl n -|_l<5n+l — (2/1 ^ra)(2/2 ^ra) ®n[(2/l 

Let 8 n = 0 n /0 n -i (yes, this is the 0 n 

®n-l^n "“ and Q n —iU n — Xrj,, 

^n+i^n+i ■ (l/i b n ^Q n a n y n , 

®n+iX^+i = ( 2/2 — b n )0 n a n T n , 

0n+l = (2/1 — h n ) ( 2/2 bn)Q n Un[(2/l — b n )Y n T (2/2 frn)Xn] “I” d n ®n— 1 • 

Remark that a n 0 n _iT n and a n 0 n _iZ n are independent of n. 

These recurrence relations for Y n , and 0 n are exactly the recurrence 
relations satisfied by products of solutions of (2) at y\ and y 2 ] . Indeed, if 
®n+lfn+l (2/1 ^n)£n ®nCn—1) and (//2 b n ) Tj n (l rl f) n . \ . 

one finds for [£ n ?7n-i> Zn-iVn, ZnVn] exactly the recurrence for [T n ,y n ,0 n ] (Ac¬ 
tually, this is a recurrence of fourth order, one should work with vectors 

\f,nPn— 15 £n— lPm CnPrn £n— iVn— l] arid [Y n , y n , 0 n , 0 n _i]). 

Now, any solution of (2) is a combination of p n {x ) and g n (x) defined by 
q n (x) = f p n (t){x - t)~ 1 dp(t) = l/{^ n x n+1 ) -|-, 

■J Supp. fj, 

a useful well-known identity is 


bp, ) X n T (l)2 &n)C/n]/T ®n /■ 

which will appear in (9)), then, with 


PnQn— 1 Pn—lQn — 'Yn/'Yn—l — l/®nj 


so, f n is some combination of p n (yi) and q n ( 2 / 1 ), Vn is some combination of p n (y 2 ) and 
q n (y 2 ), and 


T 

n 

= ( 2/1 )Pn- 1 ( 2 / 2 ) 

+ 

PPn (2/1 )tfn-. l(2/2) 

+ 

79n (2/1 )Pn-1(2/2) 

+ 

(2/1 )?n-1(2/2 

Xn ~ 

= ap n — 1 ( 2/1 )p n ( 2 / 2 ) 

+ 

(dpn-i(yi)qn(y 2 ) 

+ 

7«n-l (2/1 )Pn (2/2) 

+ 

8 qn— 1 (2/1 )?n (2/2 

©n = 

= Q:Pn (2/1 )Pn (2/2) 

+ 

fiPn (//l) Qn (2/2 ) 

+ 

72?ro (2/l )Pn (2/2 ) 

+ 

8 q n (yi)qn{y 2 ) 


As we want the degrees of the left sides to remain bounded when n increases, we must 
have a = 0. Let us show that <5 = 0 as well: from p n (yi) = X n p n _i(yi) + Y n p n -i ( 2 / 2 ) , 
<C*n-lPn{yi) - @n-lX n p n -i(yi) = 0„uHPn(2/l) ~ ^ nPn-lM = [fiPn- 1 (2/1 )Qn- 1 (2/2) + 
79n-1 (2/1 )Pn-1(2/2) + 8q n -1 (yi)qn -1 (2/2)]Pn(2/1)“ [/?Pn(2/l)?n« 1 ( 2 / 2 ) + 79n(//l )Pn- 1 ( 2 / 2 ) + 
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Sq n (yi)Qn-i{y 2 )]p n -i(yi) = iPn-i(y 2 )[q n -i(yi)Pn(yi) - q n (yi)Pn-i{yi)]+ 
Sq n -l(y 2 )[qn^l(yi)Pn(yi)-qn(yi)Pn-l(yi)\ = b/Pn-lM+Sqn-lM]/= ©n-l^nTn -1 ( 2 / 2 ) , 
or 

7 P n _i(y 2 ) +Sq n -i(y 2 ) = a n Q n - 1 Y n p n _ 1 (y 2 ) = constant p n _ 1 {y 2 ) , possible only if the 
constant (with respect to n) a n Q n -\Y n = 7 , and <5 = 0. 

One finds similarly a n Q n -iZ n = /3. 

Finally, at n = 0, 0 O = Ppo(yi)qo(y 2 ) + WoMpoM = [PS(y 2 ) +lS(y 1 )]/p 0 , 
and this relation between S(y 1 ) and S(y 2 ) is exactly (8)! One has U = @ 0 , W = 

( 2/2 - ~ 7)/(2/i 0 ), and V = -{(5 + 7 )/( 2 /x 0 ). 

□ 

4. Semi-classical measures on snuls. 

Let us consider meromorphic Stieltjes functions S, corresponding therefore to dis¬ 
crete (atomic [8]) measures. From (8), or from the equivalent form / 3S{y 2 ) + 7 <S , (yi) = 

Ho U, we have a recurrence 

(3(x k ,y k )S(y k + 1 ) +7 (%k,yk)S(y k ) = H 0 U(x k ), 

showing that poles occur at some available lattice ..., y k , y k + 1 , • • • with residues (masses 
of the measure) satisfying 

P(xk, yk)n(yk+ 1 ) + i(x k , y k )p(yk) = 0, (12) 

a Pearson-like equation like this is discussed in [29,30]. 

The masses will usually not make an infinite convergent sequence, so that they 
must remain in finite number if one wants a discrete measure. This is possible only 
if some value of (3{x,y) on a lattice, say /3(xo,yo) vanishes, so that we may start the 
nonzero masses at y 1 (and put n(yo) = 0 )? an d also if some value of 7 (x,y) vanishes 
on the same lattice, say 7 (xn, y n) = 0 , so that we stop the nonzero masses at y n (and 
put n(yN+ 1 ) = 0), see [29] (the “uninteresting case”, p.655). 

For instance, if V = 0 in (8), we simply have equal masses: /3(x,y) = —~f(x,y) = 

W(x), and W must vanish at two points of the x— lattice at least. If this does not 
happen, the measure is approximated by a discrete measure with many small equal 
masses at more and more lattice points, and tends towards the limit distribution of 
these lattice points. Take for instance ( 6 ) with x c = y c = 0, Qy/A = (y/C = 1/2, and 
q = exp(?'$), then, x k = cos k6, y k = cos (k — 1/2)0, the lattice VI of [25], distributed (if 
8/tv is not rational) like dp(x) = ( 1 — x 2 )~ 1 ^ 2 dx: we recover the Chebyshev polynomials! 

This is not surprising, as these polynomials have interesting properties with respect to 
the V —operator [ 6 , 12]. 
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5. Dense discrete measures and combinations of Chebyshev polynomi¬ 
als. 


We keep Xk = cos k9, yk = cos (k — 1/2)9 , the lattice VI of [25] with 6/n not 
rational, and try to find a discrete measure with jumps at each of these yk- “Strange” 
supports, er, carriers , have been well worked [11,20,23,31,32,33], there is a reference to 
Stieltjes himself in p.202 of [18]. 

As semi-classical orthogonal polynomials do not seem to be related to dense dis¬ 
crete measures, one tries combinations of the simplest such items, i.e., Chebyshev 
polynomials: 

Theorem. Let the measure y be discrete with jumps y(yk) = 1 /(fc — 1/2) 2 at yk = 
cos((k — 1/2)9), k = ... , —2, —1, 0,1,2,..., with 9/ 7r irrational. Then, the orthonormal 
polynomials are 

Po = 


= (f2n-l + L2n-l)T n ~ (- 1 ) L^-V/CVJ , 2n _ 1 T| n _g 3 _ 1 1 + (- 1 ) L^n-i. 9/(2*)} 

\/2Tl 2 £2nP2n{£2n-l + t2n-l) 

n > 1 (13) 

where is the value of p which minimizes p9/ (27r) — [p9/(2ir)\ on p = 1,2,... ,j, 
rjj is the value of p which minimizes [p#/( 27 r)J + 1 — p9/(27 t) on p = 1,2,... ,j, 

( where [x\ is the largest integer smaller or equal than x), £j = / 3 9/[2t r) — [tjj9/ (27r)J, 
H = lVjO/(2n)\ + 1 - Vj0/(27r). 


Remark that £j and tj are positive and decreasing with j. These ^’s and jfs are 
known as denominators of remarkable rational approximants to the irrational number 
9/{ 27t) (“Nebennaherungsbriiche” in [26] § 16) and are linked to the continued fraction 
expansion of 9/(2tt). Each new § or rj is the sum of the two last ones: if £j > ij, t — 
& = ij + Vj an( i ££ = £ j ~ tj , i£ = ij ; if £j < ij, t = rji = tjj + rjj and i£ = ij — £ 3 , £i = £ 3 . 
For instance, with the golden ratio 9/(2n) = (5 1 / 2 + l)/2 = 1.6180339887498948481..., 
one encounters the... Fibonacci numbers: 


£ or Tj 


1 

2 

3 

5 


g = j9/ (27t) - \Jj9/ (27r)_ 
0.6180339887498948481 
0.2360679774997896963 

0.0901699437494742407 


= \j]9/(2n)\ + 1 - ij9/ (27t) 


0.3819660112501051518 


0.1458980337503154555 


0.0557280900008412147 


13 0.0344418537486330259 
21 

34 0.0131556174964248372 
55 


0.0212862362522081887 

0.0081306187557833515 
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With 0/(2 tt) = y/2, one finds new values of £ at 1, 3, 5, 17, 29,..., new r\ values 
at 1, 2, 7, 12, 41,... 

Remark also that the present discrete measure is definitely not semi-classical, as 
p{yk+i) / piyjk) = ((k — 1/2) /(k + 1/2)) 2 is not a rational function of Xk = cosk9 and 
Dk = cos ((k — 1/2)9), as it should have been according to (12). 

Proof of the theorem. In order to show that a form like p n = A n T n + B n T n _£ 2ri _ 1 + 
C n T n _ V2n _ 1 is valid (using Tj p | = T p ), one must show that the scalar product 


(Pm I'm) 


oo 

'y ^ y(yk)Pn(yk)Tm(yk) 

k= — o o 


1 

2 


oo 

y ^ P{tJk) \An(T n -\-m{.yk)ATn—m{.yk))~\~B n (T n J rrn —^ 2 n _ 1 (l)k)~\~T n — r n —^2 

k=—o o 


-\~C n (Tn-\-m—rj2ri_i(.yk) T Pn—m — p2r. 


Avk))_ 


vanishes for m — 0,1,..., 2n — 1, and has the value 2 n 1 /y n = 1/A n when m = 2 n. 
Let 


OO 

t p= ZZ ^yk) T p(yk) 

k= — oo 
oo 

= ZZ (k ~ 1/2) -2 cos(p(k — 1/2)9) 

k= — o o 

= 8[cos(p0/2) + cos(3p0/2)/9 + cos(5p0/2)/25 + • • •] 

= 2tt 2 (-1)^/( 2 ^)J ([ P 9/{2tt)\ + l/2-p9/(2n)) (14') 

= -27T 2 (-l) rp0/(27r)1 (\p9/( 2tt)1 - 1/2 — p9/(2tt)) , (14") 


from elementary Fourier series, where [~af| is the smallest integer larger or equal than 
an The form of (14") will sometimes be more convenient than (14'). The importance 
of using Chebyshev moments of measures with discrete masses has been shown by 
Prevost ([27,28]), remark in particular that r p does not tend to zero when p —> oo, 
as it should with absolutely continuous measures. So, we have ( p n ,T m ) = ( A n r n+m + 

■^nhi+m-^2n — l ACnTn+m — r}2n — l An'kn — mB B n T n — m — £,2n_ l m—r)2n — l )/^‘ Remark 


that To = 7r 2 , so that po = (to) 1//2 = 7r 1 . For n > 1, let us show that it is possible 

+ C n r N - V2n _ 1 = 0 for N = 


to find A n , xv n , 
1, 2,..., 2n — 1: 


B n , and C n such that A u tn + B n TN-^ 2r 


let p = 9/ (2-7t). As 0 < £ 2 n-ip ~ L&n-ipJ < Np- [Np\, (N - £ 2 n-i)p - 1 < L^pJ - 
&n-lP < L^pj - L6n-ipj < (N-&n-l)p, One has [WpJ - L&n-lpJ = [(N-^ 2 n-l)p\- 
As |_?/2n—lPj + 1 - V2n-lP < Np\ + 1 - Np < 1, (N - T]2n-l)p < [Np\ ~ |_P2n-lPj < 
-L ? /2n-ipJ + Np < 1 + (N - rj2n-i)p, one has [Np\ - |_P2n-ipJ = \(N - y 2 n-i)p\- 


-i M) 
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In summary: 

[Np\ - L&n-lpJ = l(N ~ &n-l)p\ , [Np\ ~ [V2n-lp\ = \(N ~ T ) 2n -1 )p \, (15) 

for N = 1,2,... ,2n — 1. 

So, t N -^ 2ji _ 1 = 27 r 2 (-l)L^J-L^- 1 pj (L7Vp j _ l& n _ lP \ + 1/2 - Np H 2 n-ip) = 

(_!)Le2n-ipJ Tjv _|_ 2 7 r 2 (-l)L Ar ^-^-ipJ £2n _ l5 

T N- mn -i = -27T 2 (-l)L^J-L^-ipj ( L7Vp j _ [ V2n _ lP \ - 1/2 — N P +i l2n - lP ) = 

-(-l)L 7 ? 2 —ipJtat + 27r 2 (-l)L^J-'L€ a -iPJ t2n _ 1> 

A n TN + B n TN -£ 2n _i + C n TN- V2n -i = [^n + ( — 1) L&n-ipJ — ( — 1) lV2n-i pj rAr q_ 

27t 2 (—1)L^pJ [(_ l)\-^ 2n - 1 Pie 2n -iB n + (—l)L r ? 2 n- i pJ t 2 n _iC' n ], which vanishes indeed for 
N = n±m = 1, 2,... ,2n—1 if-B n = -K n (-l)^i P } i2n _ u Cn = K n (-l)^-^e 2 n-i, 
and A n = K n (e 2 n —i + *- 277 - 1 )- 

We now have to look at (p n ,T n ) = ( A n r 2n + B n r 2n —^ 2n _ 1 + C n T 2n — r}2n _ 1 + A u tq + 

^" T 6n-l )/ z = ^-/A n . 

1. If £277 = £ 277-1 < 2n and ?? 2 n = V 2 n— i < 2n, nothing changes in the 

evaluation of A u tn + B n TN-£ 2n _ 1 + C n TN- ri2n _ 1 when we replace TV by 2n, so 
A n T 2n + B n r 2n —£ 2n _ 1 + C n r 2n - mn _ 1 = 0, and we have only to look at A n r 0 + 
B n n 2n _ x + C n r mn _ x . One has r 0 = tt 2 , r^ 2n _ x = 2tt 2 (-1) l ^- 1 ^(1/2 - e 2n -i), 
T m n -i = — 2 tt 2 (— 1 )L^ 2 r.-ipJ (1/2 - fc 2 n-i), yielding A n r 0 + B n r^ n _ x + = 

4K n ir 2 £ 2n _iL 2n _i, which must be equal to 2/A n = 2/(K n (e 2n _i + i 2n -i)), whence 
K n = [(e 2 n—l + t2n-i)/(27r 2 £2n-it2n-i)] 1 ^ 2 , and this gives (13), as we still have 
£ 2 n = ff 2 n-l and i 2n = l2n-l' 

2. £ 2 n — 2« or r/ 2n = 2n, which happens only if £ 2n _i + V 2 n-i = 2n. Then, A n To + 
B n T£ 2n _ 1 + C n T 1l2ri l + A n T 2n + B n T 2n -£ 2rl l + C n T 2n - r)2rl _ 1 = A n (ro + r 2n) + {B n + 

C , n )(r ^ 2n -1 + r ?72n-l)- 

An interesting consequence of (15) is that [£ 277 -ipJ and |_ r / 2 n-ipJ have now the 
same evenness: if we subtract the two equations of (15) with N = n = (£ 271-1 + 
? 72 n— i)/ 2, one finds -|_£277-lpJ + L r /2n-ipj = L(P2n-l - £ 2? 7 -i)p/2J - [(£ 277-1 - 
P2n-i)p/2], which is an even integer, as |+J = — |~— x] ([15] § 1.2.4). 

So, let cr = (-l)Lftn-ipJ = (_i)b72„-ipj. One has B n + C n = <rK n (e 2 n _i - i2n- 1 ), 
k 2 n — 1 + 7-Tjan-i = 27T 2 cr(l/2 - £ 2n — 1 ) - 27T 2 cr(l/2 - t 2n -l) = 27T 2 cr(l 2 n-l “ £2n-l)- 
2a. If £ 2,7 = 2 n, £ 2 n = £ 277 - 1 -*- 277 - 1 , 2np-|_2npJ = £277-ip-L£2n-ipJ-(L ? /277-ipJ + 

1 - P277-1P), SO, [2npJ = L£277-lPj + L r /2n-ipj - 1, T 2 ?7 = 27T 2 (— 1) L 2ri Pj (1/2 - 
£2 n) = 27t 2 (1/2 — £ 277), A n {jo + T 2n ) + (B n + C n ){r^ 2n _ 1 + T rj2n _ 1 ) = 

4:7t 2 K n £ 2n L 2n — 1, whence (13), as one still has t 2n = 1277-1- 
2b. If ??277 = 2n, i 2n = t 2 n-i - £277-1, L 2 ^pJ + 1 - 2 n P = [p 2n -i P \ + 1 - 
P277-1P - (£277-IP - L£2n-ipj), SO |_2rzpJ = LP277-IPJ + LC277-IPJ, T2n = 
27T J (-l) L2n/5j ( — 1/2 + 1277 ) — — 27T" j (1/2 — 1277), A n (To+T 2n ) + (B n +C n )(T£ 2rl l + 
r V 2 n-i) = 47r 2 Ar n i 2 ,7£2n-i, whence (13), as one still has £ 2n = £ 2n -i- □ 

A numerical check has been performed with P = 0/(2tt) = (5 1 / 2 + l)/2 = 
1.6180339887498948481..., the recurrence coefficients have been computed by Gautschi’s 
sti (Stieltjes, of course!) subroutine [9]. 
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For each n, one compares the computed Ti^ n 2~ n = l/(2 n a\... a n ) with the formula 
predicted from (13), i.e., ((£ 2 n-i + ^ 2 n-i)/( 8 e 2 ?i^ 2 ?i)) 1//2 - The agreement is satisfactory, 
taking into account that all the series involving /J,(yk) = 1 /{k — 1 / 2) 2 have been trun¬ 
cated to 20000 terms, so that relative errors of about 10 -3 may be expected. 

We have longer and longer intervals where a n = 1/2 and b n = 0. On this example, 
liminf^^ooa^ > 0 , although the vanishing of this liminf could have been expected 
from pure discrete (atomic) measures [ 8 ], but other results have been published about 
singular measures [7,11,13,14,18,20,23,31,32,33]. 
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